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Nf = 24: The octagons in Fig. 1 are from a 
24 flavor staggered fermion simulation on 6^ x 4 
lattices. The data points agree with the analytical 
prediction for m > 0.25. At m = 0.25 simulations 
with Nf < 8 agree with the analytic prediction. 
The deviation here should be attributed to the 
fact that l3c = 3.90(5) is a very strong coupling 
where /3 = 6/gQ does not hold anymore. 

Nf = 17: (Plusses) These data are from runs 
using the Langevin updating algorithm on Nt = 
4 lattices. The analytic formula consistently over- 
estimates the shift in /3c. This is hard to under- 
stand given that the Nf = 24 and Nf = 8 sim- 
ulations (see below) are well represented by the 
formula. We suspect the deviation is due to inte- 
gration timestep errors in the simulation. 

Nf = 8: (Diamonds) These data are also from 
runs using the Langevin updating algorithm on 
Nt = 4 and 6 lattices. The analytic formula ac- 
curately predicts the location of the transition or 
crossover point for the larger values of the quark 
mass studied. At very small m and Nt > 8 a 
number of authors [4] have seen a transition which 
may be a bulk transition. Our analytic formula 
does not predict this transition. On the other 
hand at such small mass values Eqn. 4 is not ex- 
pected to be valid anymore. 

Nf = 4: (Bursts) These simulations do not 
show a phase transition at moderate values of the 
quark mass. At small m they show a first or- 
der transition which is believed to be associated 
with chiral restoration. The location of the tran- 
sition/crossover is well predicted by Eqn. 4 down 
to m = 0.05. For #t = 4 at m = 0.073 the first 
order chiral transition switches on [4]. Surpris- 
ingly, Eqn. 4 is still valid. One can see deviation 
from the analytic formula for m < 0.025. 

Nf = 2: (Squares) Most of the Nf = 2 simu- 
lations were performed at very light values of the 
quark mass. They do not show a phase transi- 
tion; instead, they show a smooth crossover from 
a chirally broken phase to a chirally restored one. 
Nevertheless, the location of the crossover point 
is very well tracked by the analytic formula, even 
at very light values of the quark mass. 

The agreement of the data in Fig. 1 with the 
analytic prediction, especially with smaller Nf, is 
remarkable even for masses as small as m = 0.05 



or below. The fact that the data appear to lie 
on a universal curve is a signal that the fermions 
induce an effective /3 whose strength is linear in 
Nf at fixed quark mass, down to very small mass. 

4. SUMMARY 

We demonstrated that the effects of fermions 
on the finite temperature phase transition can be 
described by an induced effective plaquette term 
for masses as low as m ~ 0.05. The induced cou- 
pling is proportional to the fiavor number and is 
independent of Nt ■ The proportionality constant 
is given by a simple 1-loop formula. From the 
point of view of lattice simulations of QCD, our 
results show that dynamical quarks must be very 
light to cause interesting effects. A finite tem- 
perature simulation at some quark mass ought 
to show an induced /3 which is not given by the 
one-loop formula, before one could claim that a 
T = simulation at the same mass would be sen- 
sitive to the effects of dynamical quarks. This is 
just barely the case in contemporary dynamical 
fermion simulations. 
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When is the above approach valid? The ques- 
tion is two-fold: 1) Can the non-local effective 
action indeed be replaced by a single plaquette 
term and 2) how well does Eqn. 4 predict the co- 
efficient of this term? It is possible to have a pure 
gauge effective action in a region where Eqn. 4 is 
no longer valid. 

3. VALIDITY OF THE EFFECTIVE 
MODEL 

We studied the dependence of the finite tem- 
perature confinement-deconfinement phase tran- 
sition on Nt , Nf and m. 

The quenched phase transition at Nt = 4 is 
at l3f = 5.69(1) [4]. Introducing Nf fiavors of 
fermions with mass m will shift the transition to 
/3f ^ = 13^ - A/3. If m is such that the fermionic 
action can be considered pure gluonic at low en- 
ergies then Al3(Nf , m) = Nf Al3i(m). If, in addi- 
tion, m is large enough that the perturbative for- 
mula is valid, A/3i(m) is given by Eqn. 4. Thus 
we expect the following behavior for the shift 
Al3(Nf , m): For m >> A where the fermion and 
gluon mass scales are well separated we expect 
to see universal behavior A/S/Nf = f(m) where 
/(m) is given by Eqn. 4. For smaller m we ex- 
pect Eqn. 4 to fail quantitatively. However, it 
might happen that A/S/Nf is still some univer- 
sal function of the quark mass. Finally, when the 
fermion scale is the same order as the gauge scale 
one can no longer replace the fermions by an effec- 
tive gauge action. The shift A/S/Nf would then 
be different for different Nf , Nt ■ Measuring the 
finite temperature transition for different Nf and 
m values makes it possible to distinguish the dif- 
ferent scenarios. 

The finite temperature transition is first order 
for the pure gauge theory, and is stable under the 
inclusion of heavy fermions. At sufficiently light 
quark mass the deconfinement transition line ter- 
minates. We might still be able to track the 
crossover point as a function of Nf and m, as long 
as the fermionic spectrum remains heavy com- 
pared to the low energy gluon spectrum. 

At very small or zero quark mass (depending 
on the number of light fiavors) there is a sec- 
ond transition whose behavior is thought to be 



primarily chirally-restoring. At this transition 
the role of the fermions is fundamental and one 
would not expect the decoupling of the gluonic 
and fermionic spectrum. 

In the above consideration we had to assume 
the relation /3 = 6/gl - the induced gauge cou- 
pling is expressed through the bare continuum 
coupling gQ while in a lattice simulation one uses 
the coefficient of the plaquette term j3. j3 = 6/gQ 
should hold in the continuum, large /3 limit; one 
expects to encounter deviations when the finite 
temperature transition happens in the strong cou- 
pling (small l3) region. 

In the following we collect the available numeri- 
cal data on finite temperature phase transition for 
different Nt, Nf and mass values [1]. To compare 
simulations with different Nf and Nt values we 
translate numerical data to express the shift in 
the gauge coupling caused by one of the fermions 
only, A/3i = {(]^ - Pc^ )/Nf. Here is the 
Monte Carlo quenched critical coupling and f]^^ 
is the Monte Carlo Nf fiavor critical coupling. 

The results are collected in Fig. 1 where the 
solid line corresponds to the analytical prediction, 
Eqn. 4. 
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Figure 1. The induced gauge coupling divided 
by the number of fiavors, A/S/Nf, compared with 
the curve from Eqn. (4) as a function of quark 
mass. Data are labeled with octagons for Nf = 
24, plusses for Nf = 17, diamonds for Nf = 8, 
bursts for Nf = 4, and squares for Nf = 2. 
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Heavy dynamical fermions with masses around the cut-off do not change the low energy physics apart from a 
finite renormalization of the gauge coupling. In this paper we study how light the heavy fermions have to be to 
cause more than this trivial renormalization. 



1. INTRODUCTION 

Heavy fermions in QCD studies are "real" : 
they are present either as real physical effects 
(heavy quarks) or as the consequence of the reg- 
ularization (Wilson fermions). In all cases the 
influence of heavy fermions at low energies is ex- 
pected to be no more than some induced effective 
gauge coupling. 

When can we expect that the fermions influ- 
ence the physical spectrum in a non-trivial way 
and when can we just replace them with an ef- 
fective local gauge action? The answer obviously 
depends on the physical processes we are inves- 
tigating. Heavy fermions are always present in 
the spectrum, unless their mass is above the cut- 
off, but if the low lying gauge and light quark 
hadronic spectrum is much below the energy level 
of the heavy fermions they will not directly influ- 
ence the low energy spectrum. 

This paper is the summary of Ref. [1] where 
we addressed the above questions in detail. 



2. THE INDUCED EFFECTIVE GAUGE 
ACTION 

The fermions' induced gauge coupling can be 
calculated by evaluating a 1-loop graph if the 
fermions are heavy [2] [3]. Consider the lat- 
tice regularized model of Nf fundamental (Wil- 
son) fermions interacting with SU(3) gauge fields, 
whose action is 

S' = /3VTr([/p)+ V^„A'„™[[/]^™, (1) 

n,jj, n,m 



where 

+(r + lM^^n-^,^)- (2) 

r = 1 corresponds to the usual Wilson fermion 
formulation while r = describes Nf = 16 x Nf 
staggered fermions. Integrating out the fermions 
we obtain the effective gauge action of the form 

Seff=S,-Tr\nK[U]. (3) 

Using the continuum representation of the gauge 
field Un/j. = e"^^''(") one can express 5'^^™ in 
terms of the continuum fields j4^(n) as the sum 
of one loop diagrams. The leading term of the 
effective action is the usual continuum gauge ac- 
tion \FavFav where the coefficient l/f/g can be 
calculated by evaluating two 1-loop graphs. 

The result is given by a four-dimensional lattice 
integral 

where S{p) is the lattice fermion propagator 
S~Hp) = 5^-''E^cos(p^)-iX^ 7Msin(p^) and 
QiPii) = irsm{pf,) + jf, cos(p^). 

The effective action has additional terms con- 
taining more derivatives and/or external gluon 
legs. These graphs are multiplied by nega- 
tive powers of m and are suppressed for heavy 
fermions [2]. In the limit where the higher order 
terms can be neglected the effective action is in- 
deed a pure gauge action with bare coupling con- 
stant given by Eqn. 4. In terms of the plaquette 
action lattice model it corresponds to an effective 
plaquette term with coefficient A/3 = G/gl- 



